Molecular dynamics recipes for genome research

Molecular dynamics (MD) simulation allows one to predict the time evolution of a system of interacting particles. It is widely used in physics, chemistry and biology to address specific questions about the structural properties and dynamical mechanisms of model systems. MD earned a great success in genome research, as it proved to be beneficial in sorting pathogenic from neutral genomic mutations. Considering their computational requirements, simulations are commonly performed on HPC computing devices, which are generally expensive and hard to administer. However, variables like the software tool used for modeling and simulation or the size of the molecule under investigation might make one hardware type or configuration more advantageous than another or even make the commodity hardware definitely suitable for MD studies. This work aims to shed lights on this aspect

Multifaceted enrichment analysis of RNA-RNA crosstalk reveals cooperating micro-societies in human colorectal cancer

Alterations in the balance of mRNA and microRNA (miRNA) expression profiles contribute to the onset and development of colorectal cancer. The regulatory functions of individual miRNA-gene pairs are widely acknowledged, but group effects are largely unexplored. We performed an integrative analysis of mRNA–miRNA and miRNA–miRNA interactions using high-throughput mRNA and miRNA expression profiles obtained from matched specimens of human colorectal cancer tissue and adjacent non- tumorous mucosa. This investigation resulted in a hypernetwork-based model, whose functional back- bone was fulfilled by tight micro-societies of miR- NAs. These proved to modulate several genes that are known to control a set of significantly enriched cancer-enhancer and cancer-protection biological processes, and that an array of upstream regulatory analyses demonstrated to be dependent on miR-145, a cell cycle and MAPK signalling cascade master regulator. In conclusion, we reveal miRNA-gene clusters and gene families with close functional relationships and highlight the role of miR-145 as potent upstream regulator of a complex RNA–RNA crosstalk, which mechanistically modulates several signalling path- ways and regulatory circuits that when deranged are relevant to the changes occurring in colorectal carcinogenesis

MicroRNA co-expression networks exhibit increased complexity in pancreatic ductal compared to Vater’s papilla adenocarcinoma

iRNA expression abnormalities in adenocarcinoma arising from pancreatic ductal system (PDAC) and Vater’s papilla (PVAC) could be associated with distinctive pathologic features and clinical cancer behaviours. Our previous miRNA expression profiling data on PDAC (n=9) and PVAC (n=4) were revaluated to define differences/ similarities in miRNA expression patterns. Afterwards, in order to uncover target genes and core signalling pathways regulated by specific miRNAs in these two tumour entities, miRNA interaction networks were wired for each tumour entity, and experimentally validated target genes underwent pathways enrichment analysis. One hundred and one miRNAs were altered, mainly over-expressed, in PDAC samples. Twenty-six miRNAs were deregulated in PVAC samples, where more miRNAs were down-expressed in tumours compared to normal tissues. Four miRNAs were significantly altered in both subgroups of patients, while 27 miRNAs were differentially expressed between PDAC and PVAC. Although miRNA interaction networks were more complex and dense in PDAC than in PVAC, pathways enrichment analysis uncovered a functional overlapping between PDAC and PVAC. However, shared signalling events were influenced by different miRNA and/or genes in the two tumour entities. Overall, specific miRNA expression patterns were involved in the regulation of a limited core signalling pathways in the biology landscape of PDAC and PVAC

Systematic Analysis of Mouse Genome Reveals Distinct Evolutionary and Functional Properties Among Circadian and Ultradian Genes

In living organisms, biological clocks regulate 24 h (circadian) molecular, physiological, and behavioral rhythms to maintain homeostasis and synchrony with predictable environmental changes, in particular with those induced by Earth's rotation on its axis. Harmonics of these circadian rhythms having periods of 8 and 12 h (ultradian) have been documented in several species. In mouse liver, harmonics of the 24-h period of gene transcription hallmarked genes oscillating with a frequency two or three times faster than circadian periodicity. Many of these harmonic transcripts enriched pathways regulating responses to environmental stress and coinciding preferentially with subjective dawn and dusk. At this time, the evolutionary history of genes with rhythmic expression is still poorly known and the role of length-of-day changes due to Earth's rotation speed decrease over the last four billion years is totally ignored. We hypothesized that ultradian and stress anticipatory genes would be more evolutionarily conserved than circadian genes and background non-oscillating genes. To investigate this issue, we performed broad computational analyses of genes/proteins oscillating at different frequency ranges across several species and showed that ultradian genes/proteins, especially those oscillating with a 12-h periodicity, are more likely to be of ancient origin and essential in mice. In summary, our results show that genes with ultradian transcriptional patterns are more likely to be phylogenetically conserved and associated with the primeval and inevitable dawn/dusk transitions

Stepwise analysis of MIR9 loci identifies miR-9-5p to be involved in Oestrogen regulated pathways in breast cancer patients

miR-9 was initially identified as an epigenetically regulated miRNA in tumours, but inconsistent findings have been reported so far. We analysed the expression of miR-9-5p, miR-9-3p, pri-miRs and MIR9 promoters methylation status in 131 breast cancer cases and 12 normal breast tissues (NBTs). The expression of both mature miRs was increased in tumours as compared to NBTs (P < 0.001) and negatively correlated with ER protein expression (P = 0.005 and P = 0.003, for miR-9-3p and miR-9-5p respectively). In addition, miR-9-5p showed a significant negative correlation with PgR (P = 0.002). Consistently, miR-9-5p and miR-9 3p were differentially expressed in the breast cancer subgroups identified by ER and PgR expression and HER2 amplification. No significant correlation between promoter methylation and pri-miRNAs expressions was found either in tumours or in NBTs. In the Luminal breast cancer subtype the expression of miR-9-5p was associated with a worse prognosis in both univariable and multivariable analyses. Ingenuity Pathway Analysis exploring the putative interactions among miR-9-5p/miR-9-3p, ER and PgR upstream and downstream regulators suggested a regulatory loop by which miR-9-5p but not miR-9-3p is induced by steroid hormone receptor and acts within hormone-receptor regulated pathways

Supervised classification of facial expressions

Over the last decade, the statistical analysis of facial expressions has become an active research topic that finds potential applications in many areas. As the expression plays remarkable social interaction, the development of a system that accomplishes the task of automatic classification is challenging. In this work, we thus consider the problem of classifying facial expressions through shape variables represented by log-transformed Euclidean distances computed among a set of anatomical landmarks

Two single nucleotide polymorphisms in the MICA gene and sMICA plasma levels are associated with hepatocellular carcinoma development in an Italian population of HCV related liver cirrhosis

Background & Aims: We investigated the relationships between MICA polymorphisms, sMICA levels and hepatocellular carcinoma (HCC) risk in HCC patients with chronic hepatitis C virus (HCV) infection. Methods.154 HCV-related HCC cases, 93 HCV-related liver cirrhosis (LC) cases and 244 healthy controls were genotyped using KASPTM SNP method. Levels of plasma soluble MICA (sMICA) were measured in 132 HCC, 90 LC patients and in 78 controls. Results. Genotyping of MICA rs2596542 showed that G/G genotype was significantly more frequent in HCC than in controls and in HCC than in LC patients. As for MICA rs2596538 allele C and C/C genotype were significantly more frequent in HCC than in controls and in HCC than in LC cases. These results demonstrate that MICA rs2596542 and rs2596538 are associated with the risk of developing HCV-related HCC in a Sicilian population. The presence of both rs2596542 and rs2596538 among LC patients raised the risk of HCC by 3.94 times. In addition, by a machine learning classifier, we found that age coupled with either rs2596542 or rs2596538 was an important discriminating factor between LC and HCC patients. Finally, sMICA levels significantly increased during HCV-related liver disease progression, and a significant relationship with either rs2596542 and rs2596538 genotypes and sMICA plasma levels was found between LC and HCC patients. Conclusions. MICA rs2596542 and rs2596538 variants are associated with the risk of developing HCV-related HCC in an Italian population

TRIM8-driven transcriptomic profile of neural stem cells identified glioma-related nodal genes and pathways

Background We recently reported TRIM8, encoding an E3 ubiquitin ligase, as a gene aberrantly expressed in glioblastoma whose expression suppresses cell growth and induces a significant reduction of clonogenic potential in glioblastoma cell lines. Methods we provided novel insights on TRIM8 functions by profiling the transcriptome of TRIM8-expressing primary mouse embryonal neural stem cells by RNA-sequencing and bioinformatic analysis. Functional analysis including luciferase assay, western blot, PCR arrays, Real time quantitative PCR were performed to validate the transcriptomic data. Results Our study identified enriched pathways related to the neurotransmission and to the central nervous system (CNS) functions, including axonal guidance, GABA receptor, Ephrin B, synaptic long-term potentiation/depression, and glutamate receptor signalling pathways. Finally, we provided additional evidence about the existence of a functional interactive crosstalk between TRIM8 and STAT3. Conclusions Our results substantiate the role of TRIM8 in the brain functions through the dysregulation of genes involved in different CNS-related pathways, including JAK-STAT. General significance This study provides novel insights on the physiological TRIM8 function by profiling for the first time the primary Neural Stem Cell over-expressing TRIM8 by using RNA-Sequencing methodology

Histone variant macroH2A1 rewires carbohydrate and lipid metabolism of hepatocellular carcinoma cells towards cancer stem cells

Hepatocellular carcinomas (HCCs) contain a sub-population of cancer stem cells (CSCs) that are responsible for tumor relapse, metastasis, and chemoresistance. We recently showed that loss of macroH2A1, a variant of the histone H2A and an epigenetic regulator of stem-cell function, in HCC leads to CSC-like features such as resistance to chemotherapeutic agents and growth of large and relatively undifferentiated tumors in xenograft models. These HCC cells silenced for macroH2A1 also exhibited stem-like metabolic changes consistent with enhanced glycolysis. However, there is no consensus as to the metabolic characteristics of CSCs that render them adaptable to microenvironmental changes by conveniently shifting energy production source or by acquiring intermediate metabolic phenotypes. Here, we assessed long-term proliferation, energy metabolism, and central carbon metabolism in human hepatoma HepG2 cells depleted in macroH2A1. MacroH2A1-depleted HepG2 cells were insensitive to serum exhaustion and showed two distinct, but interdependent changes in glucose and lipid metabolism in CSCs: (1) massive upregulation of acetyl-coA that is transformed into enhanced lipid content and (2) increased activation of the pentose phosphate pathway, diverting glycolytic intermediates to provide precursors for nucleotide synthesis. Integration of metabolomic analyses with RNA-Seq data revealed a critical role for the Liver X Receptor pathway, whose inhibition resulted in attenuated CSCs-like features. These findings shed light on the metabolic phenotype of epigenetically modified CSC-like hepatic cells, and highlight a potential approach for selective therapeutic targeting